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ABSTRACT 

Class I methanol masers are believed to be produced in the shock-excited environment around star-forming 
regions. Many authors have argued that the appearance of various subsets of class I masers may be indicative 
of specific evolutionary stages of star formation or excitation conditions. Until recently, however, no major 
interferometer was capable of imaging the important 36 GHz transition. We report on Expanded Very Large 
Array observations of the 36 GHz methanol masers and Submillimeter Array observations of the 229 GHz 
methanol masers in DR21(OH), DR21N, and DR21W. The distribution of 36 GHz masers in the outflow of 
DR21(OH) is similar to that of the other class I methanol transitions, with numerous multitransition spatial 
overlaps. At the site of the main continuum source in DR21(OH), class I masers at 36 and 229 GHz are found 
in virtual overlap with class II 6.7 GHz masers. To the south of the outflow, the 36 GHz masers are scattered 
over a large region but usually do not appear coincident with 44 GHz masers. In DR21W we detect an "S- 
curve" signature in Stokes V that implies a large value of the magnetic field strength if interpreted as due to 
Zeeman splitting, suggesting either that class I masers may exist at higher densities than previously believed or 
that the direct Zeeman interpretation of S-curve Stokes V profiles in class I masers may be incorrect. We find 
a diverse variety of different maser phenomena in these sources, suggestive of differing physical conditions 
among them. 

Subject headings: ISM: molecules — magnetic fields — masers — radio lines: ISM — stars: formation 



1. INTRODUCTION 

Methanol masers are often found in star-forming regions. 
There are two sets of transitions seen to produce methanol 
masers. Class I methanol masers (most importantly the 36 
and 44 GHz transitions) are believed to be collisionally ex- 
cited, while class II masers (incl uding the 6.7 and 12 GHz 
transitions) are radiatively excited (ICragg et alJll992h . Class I 
and class II methanol masers are so metimes both fou nd in as- 
sociation with the same source (e.g., Slvsh et al.l 19941) . but the 
two classes of masers are very rarely seen at the same velocity 
or in close (subarcsecond) spatial overlap. 

Class I methanol masers, in which shocks dominate over in- 
frared radiation, have often been assumed to be tracing an ear- 
lier evolutionary state of st ar formation than class II methano l, 
water, or OH masers (e.g.,|Ellingsen 20061: iBreen et al.|[2010t) . 
Subcategorization of class I mas ers by physical conditions 
may be possible ( ISobolevI 11993b . leading some authors to 
speculate that line intensity ratios amon g the class I maser s 
may be a proxy for evolutionary stage dPratap et al.l 120081) . 
However, class I maser studies have traditionally been biased 
towards regions hosting other tracers of star formation, and 
the cluster environments in which class I masers are found are 
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usually quite complex, calling into question traditional mod- 
els ofjheerolu^nar^timeline of class I masers (Section 4.4 
of lVoronkov et alj|2010bl and references therein). 

Furthermore, class I masers have typically been observed 
with single-dish telescopes, which can identify whether or 
not a particular class I transition produces masers in a re- 
gion (and how bright they are) but do not have the resolu- 
tion to determine their location relative to masers in other 
transitions. Given the complex environments associated with 
clustered star formation, high angular resolution is required to 
ident ify the relations be tween masers and excitation sources 
(e.g., I Amy a et all 120091) and be tween multiple transitions of 
methanol d Voronkov et al. 2006). Higher angular resolution is 
also necessary to understand the physical conditions that pro- 
duce masers in e ach of the class I transitions, which may not 
be ide ntical (e.g.. lMenteiJl99U|Johnston et alJ1992tlSoboievl 
[19931) . 

These concerns motivated Pratap et al. (2008) to do an 
unbiased single-dish search for class I methanol masers 
in nearby molecular clouds, resulting in the detection of 
new class I maser features. Several sites within these 
clouds host previously known 44 GHz methanol masers, 
many of which have been mapped interferometrically (e.g., 
Mehringer & Mentenll997l:lKogan & Slvshll998HKurtz et alj 
20041) . Imaging the 36 GHz mase rs, the other bright tran- 
sition seen in numerous sources dHaschick & Baanl 119891; 
Berulis et all [l990l iLiechti & Wilson! 119961: iPratapetall 
20081) . has heretofore not been possible due to the lack of 
interferometers operating at this frequency. However, re- 
cent upgrades to the Australia Telescope Compact Array and 
the Expanded Very Large Array (EVLA) are allowing the 
first arcse cond-resolution images of 36 GHz masers to be 
produced dSarma & Momiiadl2009T: iSiouwerman et al.ll2010b 
iVoronkov et alj|2010al) . In this Letter, we report on the first 
EVLA maps of the 36 GHz masers in the DR21 star-forming 
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complex. 

2. OBSERVATIONS AND DATA ANALYSIS 

The EVLA was used to observe the 36.169 GHz 4_i —¥ 3 E 
line of methanol in DR21(OH), DR21W, and DR21N on 2010 
May 26. The array consisted of the 20 telescopes outfitted 
with Ka-band receivers. The EVLA was in its most compact 
(D) configuration, providing a synthesized beamwidth of ap- 
proximately 2."1 x 1"4. All three sources were observed in 
dual circular polarization centered on a fixed sky frequency of 
36.1731 GHz and correlated with the new WIDAR correlator. 
The 4 MHz observing bandwidth was divided into 256 spec- 
tral channels, giving a velocity coverage of 33 kms" 1 with 
a channel spacing of 0.13 kms -1 . Conversion from sky fre- 
quency to LSR velocity was performed with the assistance of 
the EVLA Online Dopset toofl Total on-source observing 
time was ~ 12 min per source. Typical single-channel noise 
levels were < 20 mJy beam" 1 near the center of the field 

Data reduction was carried out using the NRAO Astro- 
nomical Image Processing System (AIPS). The flux scale 
was set using 3C48. Complex gain calibration was done us- 
ing J2048+4310. EVLA transition issues precluded accurate 
bandpass calibration, although it could be determined from 
calibrator data that baseline phases were flat across the entire 
observing band and amplitudes were flat across about 75% of 
the observing band. Thus, our main results, including esti- 
mates of the positions of the detected masers, should not be 
significantly affected, although bandpass effects may lead to 
amplitude errors of a few percent over most of the observed 
velocity range. Data were imaged, and flux densities were 
corrected for primary beam attenuation. 

Minor velocity errors may be introduced by two effects. 
First, the uncertainty in the rest frequency of t he observed line 
(36 169.265 ±0.030 MHz: lMiiller et al.ll2004l) is equivalent to 
a velocity uncertainty of 0.25 kms -1 (2 channels). Second, 
observations were performed at a fixed sky frequency, since 
Doppler tracking was unavailable. The sky frequency asso- 
ciated with the rest frequency drifted by about 0.015 MHz 
over the course of the observations, possibly resulting in an 
insignificant spectral broadening of maser features. Since 
the net uncertainty in determining the central velocity due 
to fixed-frequency observing (< a fewxO.01 kms" 1 ) is only 
a fraction of the channel separation (and much smaller than 
a typical maser linewidth), the data were not corrected with 
CVEL. 

We present 229.758 GHz 8_i — » 7 E maser and 226 GHz 
continuum SMA data for all of our sources (P. Pratap et al., 
in preparation). The angular resolution of these observations 
is - l."2 x 0"9. The bandwidth of the 229 GHz observa- 
tions covered the velocity range from -50 to +8 1 km s -1 with 
a spectral resolution of 0.13 kms -1 . The flux densities of 
the 229 GHz masers were corrected for primary beam atten- 
uation. The rms noise near field center was approximately 
0.15 Jy beam" 1 for DR21(OH) and DR21N and 0.3 Jy beam" 1 
forDR21W. 

When possible, we compare our data against three other 
class I methanol maser transitions from the literature: 44 GHz 
(7 6i A + ), 84 GHz (5_j -> 4 E), and 95 GHz (8 -> 7i A + ). 
For DR21N, we also used archival 44 GHz data taken with the 
VLA in B -configuration on 2006 Jun 16. 

9 |http : / /www . via ■ nrao . edu/ astro/guides /dopset / 



3. RESULTS 

Parameters of detected masers are listed in Table Q] We de- 
tect 49 36 GHz maser features in the three observed DR21 
sources: 21 in DR21N, 23 in DR21(OH), and 5 in DR21W. 
Figures Q] and [2] show the locations of the 36 GHz masers 
relative to other class I methanol transitions mapped with 
interferometric resolution. The LSR velocities of detected 
masers are consist ent with previous single-dish observations 
dPratap et al.ll2008l) . 

Spectra of the integrated flux of individual 36 GHz maser 
features are often not well fit by a single Gaussian. In some 
cases multiple spectral peaks are evident, while in others 
there is evident skewness in the spectral peak. Both of these 
cases imply the existence of structure on scales smaller than 
our synthesized beam size (1-2"). Similar effects are seen 
in the 44 GHz transition, in which iKurtz et all (120041) de- 
tect 17 masers in DR21(OH) (several with multiple spectral 
peaks) using 10 m inutes of VLA D-configuration data, while 
lAravaetal.1 ( 120091) detect 49 masers using a much longer inte- 
gration of C-configuration data. When individual maser spec- 
tra in our data do show a single clear peak, we find that typi- 
cal linewidths are on the order of 0.15 to 0.35 kms -1 . There 
is also a signature of very weak, extended methanol emission 
in some sources, but the lack of proper bandpass calibration 
tempers our confidence in this result. 

There are 14 arcsecond-scale overlaps of a 36 GHz and a 
229 GHz maser at approximately the same velocity ( | V36 — 
V229I ^ 0.4 kms" 1 ). Of these 14 associations, the 36 GHz 
maser is brighter in 7, and the 229 GHz maser is brighter in the 
other 7, indicating that either of the two 7_i — > (J- l)oE lines 
may be brighter (with the caveat that the 36 and 229 GHz ob- 
servations did not occur simultaneously and that maser vari- 
ability is possible). A simi lar phenomenon is see n in the 
Jo (J- l)i A + series, where iKalenskii et al.1 (119941) find that 
while the integrated 95 GHz flux is in general higher than that 
at 44 GHz, the integrated fluxes of sources whose spectra con- 
sist of narrow-velocity 44 and 95 GHz features (where the two 
transitions presumably produce single masers in close spatial 
association) are comparable in the two transitions. Further 
high-resolution data will be necessary to examine the relation 
between the flux densities of different maser transitions seen 
in close spatial association with each other. 

3.1. DR21(OH) 

The most striking feature of the masers in DR21(OH) is 
the central region of the source, where numerous masers in 
all mapped transitio ns occur in an appr oximately elliptical 
structure (Figure Q). lArava et al.l (120091) identify more than 
30 masers at 44 GHz associated with two shocks in a bipolar 
outflow. To within the positional uncertainty of our obser- 
vations, all 14 of the 3 6 GHz masers in th e outflow region 
are coincident with the lAraya etaD (120091) 44 GHz masers. 
The 84 and 95 GHz masers are all located within the out- 
flow region, although th e lower angular resolu tion of the 
obse rvations (8."3 x 6."2. i Batrla & Mentenl ll988l and 6."5 x 
4"4. lPlambeck & Mentenl 19901 respectively) precludes more 
detailed analysis of multitransition positional coincidences. 
Many of the 229 GHz masers align with 36 GHz masers on 
the western side of the outflow, but they appear t o be m ore 
coincident with the 44 GHz masers of lArava et alj ( 120091) . es- 
pecially in the eastern part of the outflow. 

All transitions exhibit a redshift on the western side of the 
outflow and a blueshift on the eastern side. There is also a 
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FIG. 1 .— Class I methano l maser plot of DR21(0H) showing the 36 GHz (this work), 44 GHz <Arava et alEW% , 84 GHz jBatrla & Menter][T98H), 95 GH z 
jPlambe ck & Menten 199(j), and 229 GHz (P. Pratap et al., in preparation) class I transitions and the 6.7 GHz class II transition (Harvey-Smith et al] |2008t) . 
Contours indicate 226 GHz continuum (P. Pratap et al., in preparation). The angular resolution of the 84 and 95 GHz data was lower than for the 36 GHz masers 
and should therefore be taken as indicative only of the general distribution and velocity structure of the emission in these transitions, which are consistent with 
nearby class I masers in other transitions. The dotted circle indicates the primary beam of the 226 GHz continuum and 229 GHz maser observations. The inset 
box shows an enlargement of the 3"25 square box near the origin, where 6.7, 36, and 229 GHz masers appear in close proximity and in the same velocity range. 
Coordinates are relative to the EVLA pointing center 20 h 39 m 00?8,+42°22'48f'0 (J2000). 



general velocity gradient along the direction of the redshifted 
outflow to the west, with LSR velocities generally increasing 
from the center outward. Nevertheless, some of the masers on 
the northern and southern edges of the outflow have velocities 
that do not fit this pattern. The bright 36 GHz masers are 
located in the western side of the outflow, with the brightest 
36 GHz maser bei ng located very near the brightest 44 GHz 
(lAraya et al.ll2009l) and 229 GHz masers. 

Within the outflow, each detected 36 GHz maser is found 
in close spatial association with a 44 GHz maser at approxi- 
mately the same velocity. The reverse is not true; that is, there 
are numerous 44 GHz features with no corresponding 36 GHz 



detection. On the eastern side of the outflow, 36 GHz masers 
are seen only at the eastern edge, while other transitions are 
seen throughout. In the central and western portions of the 
outflow, many of the 44 GHz masers, especially those that are 
highly redshifted or blueshifted compared to the velocity of 
the largest portion of the masers in the western outflow, have 
no accompanying 36 GHz features. 

Class II masers in the 6.7 GHz transition (5i — > 6o A + ) 
only appear projected at op the bright continuum fe ature near 
the origin in Figure Q] ( Harvev-Smith e t al. 2008). Several 
229 GHz class I masers also appear, coincident with the 
6.7 GHz masers in some cases at the 0."1 level (one-tenth of 
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FIG 2.— Class I methanol maser plots of DR21W (left) and DR21N (right). Symbols are as in Figure[T] The 44 GHz masers in DR21W are from lKogan & Slvshl 
II 998ft . The 44 GHz data in DR21N were taken with the VLA in B-configuration. No 226 GHz continuum is detected in DR21W. Every 229 GHz maser in 
DR21N is found in close proximity to a 36 GHz maser. Coordinates are relative to the EVLA pointing center 20 h 38 m 55?0, +42°19'22;'0 for DR21W and 
20 h 39 m 02?0, +42°25'43"0 (J2000) for DR21N. 



the 229 GHz beam size) but offset slightly to the east and 
south. We also detect weak 36 GHz maser emission in the re- 
gion, with spots ranging from 100 to 150 mJy detected in eight 
spectral channels between -3.8 and -2.6 km s -1 , likely due to 
spectral blending of maser features that are not fully resolved 
at our resolution. We report the position of the strongest fea- 
ture in Table Q] This feature is offset even farther to the east 
and south than the 229 GHz masers, with the caveat that it is 
difficult to determine the position of the 36 GHz maser accu- 
rately given the weakness of the maser and the 2" beamsize 
of the 36 GHz observations. 

Several class I masers appear to the south (36 and 44 GHz) 
and north (44 GHz only) of the outflow. Unlike in the out- 
flow, the 36 GHz masers in the south are usually not ac- 
companied by cospatial 44 GHz masers. The northern and 
southern masers are generally redshifted and blueshifted, re- 
spectively, compared to the bulk of the masers in the outflow, 
although not all masers conform to this pattern. The south- 
ern masers are located closer to the sources DR21(OH)W, 
DR21 (OH)S, and the ridge of amm onia emission connecting 
them dMangum et alj|l992t see also lSiouwerman et al.ll2010h 
than to the outflow in DR21(OH), and therefore are unlikely 
to be associated with the main source driving the outflow in 
DR21(OH). 

3.2. DR21W 

Despite the detection of only 5 maser features, DR21W 
contains the brightest (> 80 Jy) 36 GHz maser in our sam- 
ple. This maser is the only one for which we detect an an- 
tisymmetric Stokes V pr ofile (Figure |3j, which will be dis- 
cussed further in Section |4~T1 The bright masers near the ori- 
gin in Figure [2] are coinc i dent w ith the 44 GHz masers de- 
tected by iKogan & Slvshl d 19981) to better than a synthesized 
beamwidth. The velocities of the 36 GHz masers near the ori- 



gin agree with the velocities of the brightest nearby 44 GHz 
masers to within 0.2 km s -1 or better. 

3.3. DR21N 

DR21N was firs t identified as a 36 GHz maser source by 
iPratap et"afl (12008b . We detect numerous 36 GHz masers in 
this source, mostly in the range from -6 to -4 km s -1 , on the 
eastern side of the distribution of maser emission (Figure 0. 
Several redshifted masers in the range -1 to +5 kms -1 are 
found on the western side, along with one highly blueshifted 
maser (-14 kms -1 ). All detected 229 GHz masers are found 
in close proximity to (< 0."5 from) a 36 GHz maser. In con- 
trast, no 44 GHz maser is seen within 1" of a 36 GHz maser. 

4. DISCUSSION 
4. 1 . Circular Polarization 

As shown in Figure [3] we detect an antisymmetric "S- 
curve" Stokes V signature in the spectrum of the brightest 
maser in DR21W. If interpreted as due to Zeeman splitting, 
the implied line-of-sight magnetic field strength is +58.1 ± 
6.2 (1 .7 Hz mG~ jz) mG, where z represents the Zeeman split- 
ting coeffici ent (assumed to be 1 .7 Hz mG -1 for the 36 GHz 
transition bv lSarma & Momiianll2009h . The direction of this 
component of the magnetic field would be oriented away from 
the observer. The full three-dimensional field strength would 
presumably be larger, as Zeeman splitting measures only the 
line-of-sight component of the magnetic field when the split- 
ting be tween Zeeman components is much less t han the line 
width dWatson&Wvldll200Tt IVlemmingsll2008l) . Stokes V 
S-curves were not detected in any of the other maser sources, 
with a 3 a upper limit of 73 (1.7 HzmG~'/z) mG in the bright- 
est remaining maser in DR21(OH), although we note that a 
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TABLE 1 
Detected Masers 


Source 


RA 
Offset 
(arcsec) 


Decl. 
Offset 
(arcsec) 


VLSR 
(kmr 1 ) 


Peak 
Flux Density 

(Jy) 


DR21N 


-1.21 


5.46 


6.07 


0.96 




3.30 


1.33 


0.76 


23.72 




3.35 


1.36 


5.42 


0.17 




5.10 


3.25 


-12.58 


0.31 




6.01 


-0.43 


3.22 


0.40 




6.19 


-0.26 


6.07 


0.19 




6.24 


-6.84 


-3.38 


1.63 




7.73 


10.83 


-1.31 


0.28 




8.56 


-1.12 


-2.61 


0.28 




9.63 


4.45 


-4.94 


0.25 




10.77 


3.29 


-5.71 


0.18 




11.01 


-3.54 


^4.42 


0.96 




11.28 


-2.57 


-4.42 


0.37 




11.57 


^.19 


^t.16 


0.65 




12.54 


8.44 


-5.71 


0.15 




12.75 


-0.09 


^t.42 


0.31 




13.07 


-4.16 


-4.03 


0.60 




13.30 


5.92 


-4.16 


0.36 




14.64 


-6.39 


-3.12 


10.08 




16.74 


-3.14 


-4.03 


0.22 




20.77 


11.47 


-5.20 


6.46 


DR21(OH) 


-25.86 


-24.96 


-1.48 


2.66 




-21.08 


-20.95 


-6.01 


0.36 




-18.88 


-30.34 


0.47 


2.67 




-16.99 


0.75 


0.47 


30.26 




-16.34 


-0.68 


-1.35 


1.26 




-14.39 


-36.24 


-5.62 


0.60 




-11.92 


1.25 


-0.05 


1.31 




-10.96 


1.77 


-1.35 


0.52 




-10.52 


-1.72 


-0.57 


1.49 




-9.93 


-2.98 


0.08 


10.39 




-7.10 


-0.57 


-1.61 


2.39 




-6.24 


-2.10 


0.08 


1.81 




-2.38 


-32.23 


-5.49 


0.82 




0.33 


-29.47 


-5.88 


0.54 




2.63 


-45.86 


-4.20 


7.61 




3.38 


-45.66 


-5.62 


0.86 




3.45 a 


0.78 


-3.68 


0.15 




9.52 


^1.73 


-6.53 


1.69 




13.40 


-7.31 


-5.88 


0.42 




13.94 


-2.89 


-3.55 


0.50 




15.48 


-3.87 


-4.71 


0.85 




15.49 


-4.07 


^.97 


0.54 




15.54 


-5.61 


-3.68 


0.46 


DR21W 


-9.57 


-18.45 


-4.90 


0.30 




-8.08 


-12.63 


-5.54 


0.43 




-7.78 


-13.23 


-2.31 


0.93 




-2.16 


0.69 


-2.44 


24.60 




-1.21 


-1.40 


-2.44 


82.29 



Note. 
+42°22'48"0 



Offsets are measured from 20 h 39"W8, 



(J2000) 



20 n 39 m 02?0, 



for DR21(OH), 
+42°25'43f'0 for DR21N, and 20 h 38 m 55?0, +42° 19'22f'0 for 
DR21W. 

a Brigh test channel of spectrally-broad weak emission; see 
Section l3~7l 

circular polarization fraction of ~ 0.2% (i.e., the V/I seen in 
DR21W) would not produce a detectable (> 3cr) Stokes V 
signal in DR21(OH) or in any of the weaker masers. 

However, there is significant uncertainty in the Zeeman 
splitting coefficient, which calls into question the direct Zee- 
man i nterpretation of methanol Stokes V S-curves. Uenl 
(119511) obtained a laboratory estimate of the Lande ^-factor 
for the 25 GHz J-i J\ series of transitions of methanol, 
but the obtained value is an average over several different 
transitions and, in any case, m ay not be appropr i ate fo r the 
36 GHz transition of methan ol. IVlemmings et all ( 120061) and 
IVlemmingsl (120081) used the Uenl d 19511) g-f actor to estimate 
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FIG. 3. — Top: Stokes I spectrum of a ~ 6" rectangular region around the 
brightest maser feature in DR21W. Bottom: Stokes V spectrum after a scaled 
version of Stokes I has been subtracted. The thick red curve shows the best- 
fit value of the scaled derivative o f Sto kes I. The detection of an "S-curve" 
in Stokes V, discussed in Section 14.11 is significant even if the spectra are 
Planning smoothed to reduce the ringing evident in Stokes I. 



the Zeeman splitting coefficie nt for the 6.7 GHz transition, 
and ISarma & Momjianl ([2009) followed their method in de- 
riving the Zeeman splitting coefficient at 36 GHz. A more 
recent look suggests that these calculations may have over- 
estimated the Zeeman splitting coefficients by an order of 
magnitude (W. H. T. Vlemmings 2010, private communica- 
tion), which would serve to increase proportionally the pu- 
tative magnetic field strengths reported fr om methanol Zee- 



man splitting in this work and others (e.g. 



yiemmingsll2008l: 
ISurcis et al.ll2009HSarma & Momiianll2009l l2010l)~ 

Even if fhe lSarma & Momjianl (12009b Zeeman splitting co- 
efficient for the 36 GHz transition is (approximately) correct, 
a line-of-sight magnetic field strength of 58 mG is uncomfort- 
ably high. Measurements of the magnetic field strength are of- 
ten used as a proxy for density. Using a scaling of \B\ oc n 0A1 
from the densities and magnetic fields of DR2 1 would give a 
number density around 10 10 cm -3 dCrutcherlll999l) . a value 
several orders of magnitude higher than tha t though to be 
suitable for 36 GHz maser excita tion ( e.g., iMentenl 1199 It 
Uohnston et aTlll992t ISobolev et al.ll2005l) . Scaling from val- 
ues typically f ound in OH masers instead (a few milliga uss, 
n w 10 6 cm" 3 : iPavlakis & Kylafislll996t iFish et al.ll2005l) re- 
sults in density estimates well in excess of 10 8 cm -3 , which 
is still substantially larger than theoretica l estimates (but see 
also Section 4.3 of IVoronkov et al. 2005). The discrepancy 
between the implied density and theoretical calculations ap- 
propriate for class I maser activity grows by two orders of 
magnitude if a smaller Zeeman splitting factor is assumed 
(W. H. T. Vlemmings 2010, private communication). 

If the Zeeman interpretation is correct, one possible expla- 
nation for the discrepancy is that the correlation between mag- 
netic field strength and density does not hold in (all) the envi- 
ronments of class I methanol masers, implying that the mag- 
netic energy density is much greater than the kinetic energy 
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density dSarma & Momiianl2009l) . Alternatively, it is possible 
that the high density implied by the large magnetic field is cor- 
rect, and the bright 36 GHz maser in DR21W is produced via 
a pumping scheme not normally assumed for class I masers. 
Typically, class I masers are thought to be produced by col- 
lisional excitation followed by spontaneous cascade down to 
lower energy levels, with sink photons escaping or being ab- 
sorbed by cold dust. This mechanism has an upper limit for 
density, above which the inverted transition is thermalized by 
collisions. One possible pumping mechanism at higher den- 
sities involves collisional excitation by a warm species and a 
collisiona l sink by a cooler species to produce the population 
inversion dStrelnitskiill 19841) . This mechanism requires a two- 
temperature mixture of particles, such as electrons (and/or 
ions) and neutrals, that may exist after the passage of a shock 
front. In principle, such a collisional-collisional pump can 
operate at arbitrarily high density, as long as the temperature 
difference between the two species is sustained. 

However, it is more probable that the Stokes V signa- 
ture we detect does not measure the magnetic field directly. 
IWiebe & Watsonl (1 1998b proposed a mechanism by which 
changes in the magnetic field orientation can convert linear 
polarization to circular polarization along the amplification 
path in a maser. The Stokes V signature produced by this 
mechanism can produce an S-curve that can mimic the effect 
of Zeeman splitting of a much larg er magnetic field. In the 
case of circumstellar SiO masers that lWiebe & Watsonl ( fl998T) 
consider, their proposed mechanism could operate with a 
magnetic field as weak as 0. 1 % of the magnetic field strength 
implied from i nterpreting the S-curve as purely due to Zeeman 
splitting. The Wiebe & Watson (1998) mechanism requires 
that the Zeeman splitting in units of frequency exceed both 
the stimulated emission rate and the decay rate appropriate 
for the masing transition, a condition likel y to be satisfied for 
class I methanol masers (see Section 4.2 of iWiesemever et akl 
12001 . 

Interestingly, DR21W (this work) and M8E 
dSarma & Momiianl 120091) . the two sources with de- 
tected Zeeman splitting in the 36 GHz l i ne, are also the 
two sources in which IWiesemever et al.l (120041) detected 
circular polarization in the 133 GHz (6-1 — s- 5o E) transition. 
Future high-sensitivity observations at 133 GHz might permit 
detection of an S-curve signature in Stokes V, thereby testing 
whether the magnetic field implied from a pure Zeeman 
interpretation is consistent with the values derived at 36 GHz. 
(Since the 36 GHz and 133 GHz methanol lines are both 
/-i — > (J-l)o E transitions, the same Lande g-factor would 
be applicable to each of them.) Careful laboratory measure- 
ments of the Zeeman splitting coefficients appropriate for 
the brightest methanol maser transitions would also be very 
helpful. 

4.2. A Class I/Class II Overlap 

Excitation models indicate that the physical conditions that 
are thought to produce inversion in class I transitions usu- 
ally lead to anti-inversion in class II transitions and vice versa 
dCragg et al.lll992l:ISlvsh et alj|2002l) . In that regard, the exis- 
tence of class I 36 and 229 GHz masers at the same velocity 
and location as the class II 6.7 GHz maser site is surprising. 
These masers are also coincident with the brightest continuum 
emission near the outflow, which may be an important clue to 
their origin. 

There is preced ent for overlaps betwee n class I and class II 
methanol masers. IVoronkov et al.l (j2005j) deduce spatial coin- 



cidences between the 6.7 GHz class II and 25 GHz J2 — > J\ E 
series of class I methanol masers in OMC-1, although the an- 
gular resolution of their interferometric observations is much 
po orer than either the 229 GH z data on whi ch we report or 
thelHarvev-SmithetalJ (120081) 6.7 GHz data. IVoronkov etail 
(12005b model the 6.7/25 GHz overlap as arising from a lower- 
temperature regime than is typically assumed for class II 
masers. This pumping model requires an intermixed environ- 
ment of gas and dust at a lower temperature (~ 60 K) than 
for traditional class II maser formation. The strong 226 GHz 
continuum emission in DR21(OH) li kely marks such an en- 
vironment that is rich in dust. The IVoronkov et al.l (120051) 
model did not investigate excitation at 229 GHz, but it does 
predict inversion at 36 GHz, which is also a J-\ — > (J- l)o E 
transition. The model predicts that 44 and 95 GHz masers 
should not appear at a low-temperature 6.7 GHz maser site, 
since the former are Jq — > (J-l)\ A + transi tions, while the 
latter is a (J— l)i — > 7n A + transition. Indeed . iKogan & Slyshl 
dl998l) . lKurtz et all d2004l) . and lArava etaTI d2009b do not re- 
port 44 GHz masers at this location. 

4.3. Multitransition Comparison of Class I Masers 

In each of the sources we observed, the 36 GHz masers di- 
vide into two sets: those that are found to be cospatial with 
other class I maser transitions, and those that are spatially iso- 
lated. We address each of these two groups in turn. 

In the outflow of DR2 1 (OH) and in the cluster in the north 
of DR2 1 W, we note that the distribution of the 36 GHz masers 
is similar to that of the 44, 84, 95, and 229 GHz masers. 
Additionally, the brightest masers in all of these transitions 
appear in the same location. In DR21(OH), the brightest 
masers in all class I transitions occur in a narrow velocity 
range (v'lsr = +0.3 to +0.5 kms" 1 ) at the western tip of the 
outflow. The next-brightest masers in each transition appear 
approximately 10" east of this position (i.e., 10" west of the 
origin in Figure [TJ. It is clear that these masers exist in a 
physical regime in which strong maser emission is produced 
in a large number of J-\ — > (J — l)o E and /(> — > (J— l)i A + 
transitions simultaneously. The existence of a nearby ener- 
getic source creating an outflow likely shocks the surrounding 
molecular material, pumping multiple transitions of methanol 
in the same location. 

On the other hand, all sources contain regions in which 
36 GHz masers are not found to be coincident with other 
transitions. For instance, numerous masers in both the 36 
and 44 GHz masers are found to the south of the outflow in 
DR21(OH). For the most part, the two transitions are not seen 
to be spatially coincident. Most of the 36 GHz masers and 
some of the 44 GHz mase rs in the south are seen ne ar the 
periphery of CS emission dPlambeck & MentenlfT990T) . The 
common feature of the environment of these masers is the 
lack of a nearby energetic source. As in the Sagittarius A 
region, it is possible that these masers trace molecular density 
clum ps that are shock-excited by core collision or compres- 
sion dSiouwerman et al]|2010l) . 

The 229 GHz masers appear to align more closely with 
the 44 GHz masers than with the 36 GHz masers in the 
DR21(OH) outflow. While 229 GHz masers are found to be 
coincident with 36 GHz masers on the western side of the out- 
flow, there is always an accompanying 44 GHz maser. In con- 
trast, several 229 GHz masers appear coincident with 44 GHz 
masers but without 36 GHz emission, especially on the east- 
ern side of the outflow. However, we note several caveats 
about drawing conclusions from these points. First, the angu- 
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lar resolution of the observations of both the 44 GHz and the 
229 GHz masers is higher than that of the 36 GHz observa- 
tions on which we report, and it is possible that more sensitive, 
higher-resolution observations of the 36 GHz masers will un- 
cover other features at the sites of the 44/229 GHz overlaps. 
Second, every detected 229 GHz maser in DR21N appears 
near a 36 GHz maser site (Figure |2). Third, at present the 
three DR21 sources on which we report are the only sources 
mapped in all three of the 36, 44, and 229 GHz transitions, 
and the field of view encompassed by the 229 GHz observa- 
tions in DR2 1 (OH) excludes nearly all of the masers south of 
the outflow as well as the eastern tip of the outflow itself. This 
constitutes a rather small set of sources from which to draw 
general conclusions about the properties of the 36, 44, and 
229 GHz masers relative to each other. More interferometric 
observations of both the 36 and 229 GHz transitions will be 
required to fully understand the phenomenology of either. 

5. CONCLUSIONS AND FUTURE WORK 

We have interferometrically imaged three sources in the 36 
and 229 GHz class I methanol maser lines. We find numerous 
masers in both transitions and a diversity of conditions among 
them. Notably, we identify the following three environments, 
each of which may correspond to a different set of excitation 
conditions. 

1 . The outflow in DR2 1 (OH) contains a large number of 
overlaps in the 36, 44, and 229 GHz transitions. Strong emis- 
sion is also seen at 84 and 95 GHz, although there is a lack 
of high-resolution observations of these lines. The brightest 
masers in all transitions appear at the same velocity and in the 
same location at the western tip of the outflow as traced by the 
masers. 

2. The bright continuum source in DR21(OH) is associ- 
ated with both the class II 6.7 GHz masers and class I 36 
and 229 GHz £-type masers, but the class I 44 GHz A + - 
type masers are conspicuously absent. This environment may 
be explai ned by the low-temper ature intermixed dust and gas 
model of lVoronkov et alj d2005[) . 

3. The 229 GHz transition produces detectable masers at 
a subset of the 36 GHz maser sites in DR21N. In both this 
source and in DR21(OH) well south of the outflow, the 36 
and 44 GHz masers have a similar large-scale distribution but 
are rarely found to produce a maser at the same site. 



In addition, the brightest maser in DR21W produces an an- 
tisymmetric Stokes V profile that implies a large magnetic 
field if interpreted as being due to Zeeman splitting. There 
is a very large uncertainty in the Zeeman splitting coefficient 
appropriate for the 36 GHz transition, but the implied density 
(assuming \B\ oc n a ,a ~ 0.5) greatly exceeds the range over 
which class I masers are thought to form. It is possible that 
the circular polarization in this maser feature is produced by a 
much smaller magnetic field whose orientation changes ove r 
the amplification path of the maser dWiebe & Watson|[T99l) . 
In any case, careful laboratory measurement of Zeeman split- 
ting coefficients appropriate for methanol maser transitions is 
warranted given the increasing number of Zeeman-like Stokes 
V signatures identified in methanol transitions over the past 
few years. 

In order to realize the goal of being able to identify the 
physical conditions in a variety of star-forming regions by the 
presence or absence of various methanol maser transitions, it 
will be necessary to understand a few sources in greater de- 
tail. Our results highlight the need for both increased theoret- 
ical effort and more sensitive observations of multiple class I 
transitions at higher angular resolution. In particular, higher- 
resolution maps of the 84 and 95 GHz transitions as well as 
higher-frequency maser lines may be both enlightening and 
timely in the advent of the ALMA era. A further survey of 
sources in the 36 and 229 GHz transitions, when combined 
with 44 GHz maps in the literature, may also be helpful in 
determining the range of possible methanol excitation condi- 
tions in nature. 
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